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ABSTRACT 

We present the first results from a new, deep (200ks) Chandra observation of the X-ray luminous galaxy 
cluster surrounding the powerful (L ~ 10 47 erg s _1 ), high-redshift (z = 1.067), compact-steep-spectrum 
radio-loud quasar 3C 186. The diffuse X-ray emission from the cluster has a roughly ellipsoidal shape 
and extends out to radii of at least ~ 60 arcsec (~ 500 kpc). The centroid of the diffuse X-ray emission is 
offset by 0.68±0.11" (~ 5.5 ±0.9 kpc) from the position of the quasar. We measure a cluster mass within 
the radius at which the mean enclosed density is 2500 times the critical density, ^2500 
The gas mass fraction within this radius is f gas 



of 1.02±g;f^ x 10 14 M 



129+ - 015 

U.1Z»_q 01 g. 



= 2831 jj kpc, 
This value is 

consistent with measurements at lower rcdshifts and implies minimal evolution in the f ga s(z) relation 
for hot, massive clusters at < z < 1.1. The measured metal abundance of 0.42j^QQ7 Solar is consistent 
with the abundance observed in other massive, high redshift clusters. The spatially-resolved temperature 
profile for the cluster shows a drop in temperature, from kT ~ 8 kcV to kT ~ 3 keV, in its central regions 
that is characteristic of cooling core clusters. This is the first spectroscopic identification of a cooling 
core cluster at z > 1. We measure cooling times for the X-ray emitting gas at radii of 50 kpc and 25 kpc 
of 1.7 ± 0.2 x 10 9 years and 7.5 ± 2.6 x 10 8 years, as well as a nominal cooling rate (in the absence 
of heating) of 400 ± 190 M Q year~ 1 within the central 100 kpc. In principle, the cooling gas can supply 
enough fuel to support the growth of the supermassive black hole and to power the luminous quasar. 
The radiative power of the quasar exceeds by a factor of 10 the kinematic power of the central radio 
source, suggesting that radiative heating may be important at intermittent intervals in cluster cores. 

Subject headings: quasars: individual (3C 186) - X-rays: galaxies: clusters 



1. INTRODUCTION 

Recent X-ray observations of nearby galaxy clusters 
show that powerful outbursts of their central cD galax- 
ies have imprinted a r ich variety of structures ont o the 
X-ray emitting gas (see iMcNamara fe Nulsenl (|2007l ) for a 
review). The amount of energy supplied into the cluster 
gas by these outbursts can p revent the clusters from cool- 
ing (jMcNamara et ai1l2005[) . The average power released 
in outbursts exceeds ~ 10 45 erg s _1 . This is equivalent 
to the typical radiative power of a quasar. However, the 
central cD galaxy is usually observed in a quiescent, non- 
luminous, state. The impact of a luminous quasar on the 
cluster gas has not been widely explored so far. 

On the other hand, it has been known for a long 
time that powerful radio- loud quasars are associated 
with rich galaxy environments (Smith fc Heckmanl 11990: 
lEllingson et~aTHl991at lYee fc Ellingsonlll993t ). and so we 
should be able to find luminous radio loud quasars in X- 
ray-bright clusters. Much remains unknown about the de- 
tails of how quasar activity is triggered, the way in which 
the accretion proceeds, the impact of the quasar on its 
galactic and cluster environment, and which physical pro- 



cesses are important as a function of redshift. X-ray clus- 
ters retain the history of the cD activity and studies of 
quasars in X-ray clusters can provide some answers to 
these questions. 

There have been a few successful searches for the diffuse 
X-ray cluster emission surroun ding powerful quasars and 
radio galaxies using ROSAT (ICrawford fc Fabian 19931: 
Worrall et al.l[l99l iHaiTet al.lll995lll997t ICrawford et al l 



19991: lHardcastle fc Wqrralll [l999t ISarazin et all 119991: 



Worrall et all l2000j) . lHall et all (1997D discussed the 



ROSAT X-ray data for five quasars located in X-ray clus- 
ters, and concluded that simple models for triggering the 
quasars (e.g., cooling flows, a low-pressure ISM, or low 
velocity dispersions) do not work, suggesting instead that 
mergers and strong interact ions are critical in del i vering 
the gas to the quasar host. lHardcastle fc Worralll (jl999) 
analyzed a large sample of powerful radio sources and 
found cases of extended X-ray emission that arc likely 
to be of cluster origin. However, the number of such 
cases was limited to eight by ROSAT's sensitivity and 
some of these were not confirmed by Chandra o bservations 
(|Crawford fc Fabianl [20031: 1 Worrall et al.ll2004D . Extended 
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X-ray emission can be associated alternatively with radio 
structures s uch as large scale jets, knots, hot spots, and 



lobes (e.g., Cclotti & Fabian 2004; Kataoka & Stawarz 
I2005t ICroston et al.l 120051: lHarris fc Krawczvnskil I2006D 
and disentangling the cluster emission from the other X- 
ray emission components has been challenging. The Chan- 
dra X-ray Observatorv (|Weisskopf et al.l 2000|) has both 
the spatial resolution and the dynamic range necessary 
to study diffuse X-ray emission in the vicinity of a strong 
point source, and allows the separation of cluster X-ray 
e mission from t h at of j ets and other structures. 

iBelsole et ail (|2007n studied a Chandra and XMM- 
Newton sample of twenty powerful radio sources at z > 0.5, 
including two core-dominated quasars, and found diffuse 
X-ray emission in 60% of the sample. The diffuse emission 
was faint, but the luminosities were consistent with non- 
quasar-host X-ray clusters at similar and lower redshifts, 
and the work found no difference between the cluster en- 
vironments of quasars and radio galaxies. 

We discovered a bright X-ray cluster in a pointed Chan- 
dra observation of the radio-loud compact steep spec- 
trum (CSS) quasar 3C 186 at the redshift of z = 1.067 
(|Siemiginowska et al.l l2005h . This observation, although 
only ~ 30 ksec long, provided an X-ray luminosity mea- 
surement, £(o.5-2fceV) = 6 x 10 44 erg s -1 , a cluster tem- 
perature (kT = 5.2+JJ keV) and a gas-mass fraction 



(/gas(^25oo) ~ 0-13 ± O. U» ) tnat were typical ot oti 
sive, relaxed clusters (jVikhlinin et al.l l2002t I AIL 
200 8j). The 3C 1 8 6 X-r ay cluster is more luminous than 



that were typical of other mas- 
en et al.l 



the IBelsole et al.l (|2007t) clusters, and provides a unique 
opportunity to study a luminous cluster associated with 
a quasar at high redshift. The powerful and luminous 
quasar, L^ Q i ~ 10 47 erg s _1 , is located well within the 
diffuse X-ray emission. 

Here we report the results from a 200 ks deep follow-up 
Chandra observation of 3C 186, made in order to study 
the quasar and its associated cluster in greater detail. 
Relatively few massive, relaxed X-ray clusters at z > 1 
are known. In addition, the 3C 186 cluster is one of very 
few X-ray clusters that is both associated with a quasar 
and bright enough for detailed study. We give details of 
the Chandra observation in Section 2, describe image and 
spectral analyzes in Section 3, and present a discussion of 
the results in Section 4. 

Throughout this paper we use c osmological parame - 
ters based on WMAP measurements (jSpergel et al.ll2003h : 
H a =71 km s~ 4 Mpc" 1 , Q M = 0.27, and ft vac = 0.73. At 
z = 1.067, 1" corresponds to —8.163 kpc. 

2. CHANDRA OBSERVATIONS 

3C 186 was observed with the Chandra ACIS-S CCD in 
December 2007. Due to scheduling constraints the obser- 
vation was split into four separate pointings that sum to 
a total exposure of 197 ks (see Table [T|). The quasar was 
placed on the back-illuminated CCD (S3) and was offset 
by -1 arcmin in Y coordinates to make sure that the clus- 
ter is not affected by a chip gap. The observation was 
made in VFAINT mode and full- window mode. Table [1] 
shows which CCD chips were active. The observation was 
not affected by solar flares, and the background was quiet 

1 http:/ /cxc. harvard.edu/chart/index. html 

2 We used Gaussian model instead of a Delta function to account for 



for the entire observation. The quasar is slightly offset in 
obsid 9408 with respect to the other three observations in 
Z coordinates. However, the four observations were per- 
formed with similar enough configurations that they could 
be merged together for the purpose of image analysis. Note 
that the exposure maps are flat on the scale size of the 
cluster in all four observations. Analysis was performed 
with the CIAO version 4 software u sing CALDB version 
4.2. All modeling w as done in Sherpa (|Freeman et al.ll2001l ; 
iRefsdal et all 120091). W e used the Cash and Cstat fitting 
statistics dCash 1119791) an d the Nelder-Mead optimization 
method (jNelder fe Mead II1965T ). We did not subtract the 
background, but included a background contribution in 
the model expressions for all of our model fitting. 

The events have been filtered to remove the VFAINT 
background events. The standard ±0.5 pixel randomiza- 
tion has also been removed before merging the four ob- 
servations. The ACIS image of the merged observation is 
shown in Figure I37T1 The bright quasar is located within 
diffuse X-ray emission that is visible on a scale exceeding 
-30" (250 kpc). 

3. DATA ANALYSIS AND RESULTS 
3.1. Image Analysis 

We performed 2D image analysis using the merged data 
to evaluate the spatial extent of the diffuse emission, and 
the quasar contribution to the diffuse emission within 
<5 arcsec of the core. Figure |3~T1 shows the merged ACIS- 
S images in three RGB colors representing different en- 
ergy ranges: red corresponds to 0.3-1.5 keV, green to 1.5- 
2.5 keV, and blue to 2.5-7 keV. We show a binned image 
with the standard ACIS-S pixel size of 0.492" and an im- 
age that is smoothed with a Gaussian kernel (a = 2.46"). 
The cluster emission is relatively smooth, elliptical, and 
elongated in NE-SW direction. We measure (using ds9) 
the largest extent of the emission to be ~ 128" along the 
major axis at PA = 43 deg, i.e. — 60" (> 500 kpc) dis- 
tance from the quasar. We selected the events within a 
0.5-7 keV energy range for the spatial modeling. 

Effects of the strong quasar emission have been taken 
into account in the analysis of the cluster data. Note that 
point sources in the field of view in Figure 13.11 indicate 
the size of the PSF for a source of typical strength. We 
adopted the observed quasar spectrum (a power law with a 
photon index of T = 1.9) as the input model for CHARlQ 
simulations of the quasar's PSF, and created a very high 
signal-to-noise PSF model image for use in our analysis. 

In our two-dimensional analysis we fitted both the clus- 
ter centroid and the quasar position to determine whether 
they are co-aligned. We tried circular models and also 
allowed the cluster gas to have an ellipticity. We used 
Sherpa and adopted a 2D GaussiarQ model for the quasar 
and a 2D Beta model for the cluster. We also included a 
constant model to account for the background. Because 
the exposure map is uniform across the cluster region we 
did not include the exposure map in this analysis, and so 
we worked with the count data and used Poisson statistics 
in fitting. However, we have confirmed that including the 
exposure map gives consistent parameters with the ones 
reported here. 

an unknown aspect "blur" in MARX. 
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Table 1 
Chandra OBSERVATIONS 



OBSID 


Exposure (ks) 


Date 


CCD 


Inner" 2.5-6.0" 


Outer 6-20" 


9407 


66.3 


2007-12-03 


ACIS-23567 


500.0±22.8 


1000.7±34.8 


9774 


75.1 


2007-12-06 


ACIS-23567 


613.0±25.3 


1110.6±36.9 


9775 


15.9 


2007-12-08 


ACIS-23567 


114.8±11.0 


241.5±17.4 


9408 


39.6 


2007-12-11 


ACIS-23567 


306.9±17.8 


627.9±27.3 



net counts within the energy range 0.3-7 keV in annuli centered on the quasar whose radii are given in arcscc. 



We used Sherpa to convolve the 2D model with the sim- 
ulated 2D PSF, and fitted the result to the data. The 
models and best- fit parameters are given in Table [2j 

The 2D fit finds a cluster core radius of r core = 3.51 ± 
0.31 arcsec (~28.6±2.5 kpc) and /3 = 0.48±0.17 for ellipti- 
cal models and r core = 3.06±0.25 arcsec (~ 25.0±2.0 kpc) 
and /3 = 0.48 ± 0.17 for the circular models. We note that 
there are systematic uncertainties present in the 2D anal- 
ysis. The circular and elliptical fits have small statistical 
errors, but the difference between the parameters obtained 
in circular and elliptical models, in particular the cluster 
core radius, could be larger than the la statistical errors 
reported in Tabled 

The best-fit locations of the quasar and the cluster 
centroids are offset by 0.68 arcsec, which is significantly 
greater than the corresponding uncertainty of 0.11 arc- 
sec. We centered the surface brightness profile on the 
peak of the X-ray emission, while the true cluster centroid 
is slightly offset from the peak. Figure [2] shows surface- 
brightness profiles centered on the quasar and obtained 
from the data and the 2D best-fit circular models. The 
quasar emission dominates over the cluster emission within 
about 1.5 arcsec of the core and it is unresolved. At 2.5" 
distance from the peak the quasar contributes only about 
10% of the total counts. The X-ray cluster emission domi- 
nates outside that region and exceeds the background level 
up to ~ 40" (~ 320 kpc) from the quasar core. The cluster 
is detected at 3a to ~ 285 kpc. 



3.2. Spectral Modeling of the Cluster Emission 

Based on our 2D image modeling results we assume that 
the X-ray emission outside the central 2.75" (~ 22 kpc) 
radius circular region is dominated by cluster emission 
(see analysis of the quasar contribution to the observed 
emission below). Note that 3C 186 is a compact ra- 
dio source with a linear size of 2", so any X-ray emis- 
sion associated with the radio source would be contained 
within this central region. We performed spectral analy- 
sis using each individual observation (obsid) and the cor- 
responding calibration files. To obtain a global temper- 
ature for the entire cluster we extracted the spectrum 
from each event file (4 obsids) assuming a large ellipti- 
cal region with the semi-minor and semi-major axis of 
20" and 30" at PA 315 degree and ignoring the inner cir- 
cle with radius of 2.75 arcsec dominated by the quasar. 
The background files were extracted from regions on the 

3 http:/ /cxc. harvard.edu/contrib/dcprojcct/ 



same CCD located outside the source region excluding all 
the detected point sources. The APEC thermal plasma 
model at rcdshift z = 1.067 was fitted to the four spec- 
tra simultaneously, giving a best global temperature of 
kT=5.58±8$ keV, a metallicity of 0.42±g$ Solar and 
a soft-band flux Fo. 5 _ 2 = 8.27 x 10" 14 erg sec -1 cm" 2 
(Fig. [3]). The best fit global cluster temperature is in 
agreement with the si ngle temperature measured in the 
first short observation ([Siemi ginows ka et al.ll2005h . 

In order to look for any temperature gradient we ex- 
tracted the spectra from each event file (4 obsid) assuming 
seven circular annuli centered on the quasar. The annuli 
cover the range from 2.75" to 30". We list the angular 
ranges spanned by the annuli in Table [3] All spectra were 
taken from individual observations to properly account for 
the instrumental effects. They were simultaneously fitted 
using all available counts for each annuli in the energy 
range 0.5-7 keV. We accounted for any background contri- 
bution using a complex empirical model (a combination of 
an 8th order polynomial and five Gaussian lines) that was 
first fitted to the ACIS-S "blank-sky" background data, 
(see Appendix Al for more details). We next fit this back- 
ground model to the background spectra from each obser- 
vation to check how well the model describes our data. In 
the simultaneous fit of the source and background model to 
the cluster spectra we varied only the background model 
normalization and kept all the other background model 
parameters frozen. We note that this is reasonable, as the 
fraction of background counts is typically lower than a few 
percent and exceeds 23% only in the outermost annulus. 

Table [3J shows the best fit parameters for the applied 
model. We assumed an APEC model, included a correc- 
tion for the Galactic absorption with an equivalent Hy- 
drogen column density of Njj = 5.64 x 10 20 cm" 2 , and 
assumed cluster metal abundances of 0.3 in respect to So- 
lar. We developed a deproject modeE in S herpa based 
on the description in iFabian et al.l (11981 ) and iKriss et al.l 
CEM). We assume a spherical geometry and the radial 
sizes of the shells given by each individual annulus. The fit 
starts at the outermost annulus and proceeds towards the 
innermost one taking into account the contributions from 
outer annuli to the fitting of the inner one. The best-fit 
model gives the deprojected temperatures, normalizations, 
and densities listed in Table El We plot deprojected tem- 
perature, density and entropy profiles in Figure 01 These 
profiles are consistent with the cluster having a cooling 
core. 
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26.0 24.0 22.0 7:44:20.0 18.0 16.0 14.0 12.0 10.0 

RA (J2000) 




26.0 24.0 22.0 7:44:20.0 18.0 16.0 14.0 12.0 10.0 

RA (J2000) 

Fig. 1. — RGB color Chandra ACIS-S images of the 3C 186 X-ray cluster. The four individual observations have been merged 
into one image. Colors represent different energies: red:0.3-1.5 keV, green:1.5-2.5 keV and blue:2.5-7 keV. Top: Images binned to 
ACIS-S pixels with the standard size of 1 pixel=0.492 arcsec. Bottom: The image has been smoothed by a Gaussian function of 
a = 2.46". 
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Table 2 

Best Fit 2D Model Parameters' 1 



Parameter Circular Models' 5 Elliptical Models b Units 



qso.fwhm 


0.43t0.02 


0.26 to. 02 


arcscc 


qso.xpos 


4072.82^0.01 


4072.82 to. 01 


physical 


qso.ypos 


3946. i4+ 0.01 


3946.43 to.Ol 


physical 


qso.cllip 




0.115 tO. 097 




qso.theta 




1.82 to. 45 


radians 


qso.ampl 


2018.5l^|;| 


1497.1 tl38.6 


counts 


clus.rO 


3. 03^0. 25 


3.58 to. 31 


arcsec 


clus.xpos 


4072.98l0.21 


4072.96 to. 21 




clus.ypos 


3945.05to.21 


3945.06 to.22 




clus.cllip 




0.300 to. 015 




clus.theta 




2.41 t0. 03 


radians 


clus.ampl 


0.42^0. 04 


0.44 to. 04 


counts 


clus.beta 


0.48 to. 17 


0.48 t0.17 




bgnd.cO 


0.0023^0.0002 


0.0024 1 0.0001 


counts 



a Fit was performed on non-background subtracted image and the background model was included as a part of the model component: qso - 2D 

gaussian model; bgnd - 2D constant background model 
clus - 2D beta model: f(x, y) = f(r) = A(l + ( ^) 2 ) -W+ - 5 ; 

r(x, y) = sjxft{\ - e) 2 + y/ 2 /(l - e); 
xl — (x — x ) cos 8 -f (y — y ) sin 8, yl — (y — y ) cos 8 — (x — x ) sin 8 8- the angle of cllipticity, e - cllipticity. 
b la uncertainties are shown for one interesting parameter. 
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Fig. 2. — Top panel: Surface brightness profiles (in units of counts per pixel) centered on the quasar extracted from the data 
and the 2D circular models (see Table 2 for the best fit parameters). The location of the quasar centroid in physical coordinates 
is marked in the upper right corner. The surface brightness profile from the ACIS-S data is shown as black points with errorbars: 
data is binned to ensure at least 200 counts per bin. The model profiles are marked with solid lines: red - full model (qso + clus+ 
bgnd as defined in Table 2) convolved with ChART PSF; green - quasar (gauss2d only) model convolved with PSF (ie no cluster 
or background component), blue - cluster (beta2d), brown - background (cont2d). Bottom panel: Residuals between the full 
model and the data in units of a. Physical pixel size is equal to 0.492", the size of the ACIS-S pixel. 
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Fig. 3. — The confidence levels for the temperature and abundance parameters obtained in a global model fit to the cluster 
described in Sec. 13.21 The location of the best fit parameters is marked by a cross. Three levels la, 2a 3a contours are shown. 



We followed iRussell et al.1 (l2010h . using the simulated 
Chandra PSF to understand the contribution from the cen- 
tral quasar to each individual annulus and the effects as- 
sociated with a possible contamination of the cluster spec- 
trum. We used CHART to simulate the quasar (point 
source) photons scattered by the Chandra mirrors, and 
MARX to project them onto the ACIS-S detector at the 
exact pointing as in the Chandra observations of 3C 186. 
In the CHART simulations we assumed the quasar photon 
flux to be described by the best model parameters fit to 
the spectrum extracted from a r = 1.75" circular region 
centered on the quasar, i.e. an absorbed power law with 
Njj and T = 1.9. The simulated quasar image reflects 
both the scattering of the Chandra mirrors and an addi- 
tional "blur" (parameter set to 0.3 in MARX) due to the 
standard dither and the aspect uncertainty. In Figure [5] 
we show a number of counts from the simulated quasar 
image in comparison to a number of counts detected by 
Chandra in each annulus assumed for the spectral model- 
ing. The fraction of the total observed counts that can be 
associated with the quasar is shown in the bottom panel 
of the figure. It exceeds 10% only in the first innermost 
annulus. 

In order to check the effects of quasar contamination on 
the spectral results, additional fits were performed wherein 
we included appropriately normalized power law model 
components to the cluster model for each annulus (sec Ta- 
ble 2]). The resulting temperatures and normalizations are 
in excellent agreement with the deprojected model results 
shown in Table 3. The only noticeable difference is for 
the innermost annulus, where the best fit temperature of 
2.54^0 57 is marginally lower (but still consistent within la- 
errors) . We conclude that the quasar contribution does not 
significantly impact the deprojected fit parameters for the 
cluster observation. 

Given the flux for the models in Table [3] and applying 
the appropriate K-correction, we find a total cluster lumi- 
nosity L .5_2 kcv = 4.6 ± 0.2 x 10 44 erg s _1 . 



3.3. Cluster Mass 



We have measured gas mass and total mass profiles for 
the cl uster using the Monte Carlo method of I Allen et al.l 
(2008]). This analysis is also used to determine gas density 
and cooling time profiles for the cluster, which are shown 
in Figure [Bj The density results are consistent with those 
obtained directly from the spectral fitting. 

Our mass analysis uses a parameterized iNavarro et al.l 
(|l995l 11997! ) (NFW) model, which is fitted directly to 
the observed cluster surface brightness profile and depro- 
jected temperatur e profile (see Eq. l and the description 
of the method in lAllen et aTll2008f ). The best-fit NFW 
model has a concentration parameter c = 7.4l?,'3, scale 
radius r s = I20I4Q kpc, and equivalent velocity disper- 
sion a = 780^gQ km s _1 , with x 2 = 7.9 for five degrees of 
freedom. For these parameters we calculate the radius at 
which the mean enclosed mass density is 2500 times the 
critical density of the Universe at the redshift of the clus- 
ter, r 2 5oo = 2831^ kpc. Thi s results agrees well with the 
results of lAllen et al.l ()2008l ) based on the earlier Chan- 
dra observation. However, the statistical uncertainties on 
?"2500 are improved by a factor of ~ 4. 

The total mass within r 2 5oo is M2500 = l-02+ ; 14 x 
10 14 M Q . (Our 68 per cent error bars on M2500 also ac- 
count for the uncertainty in r25oo-) The measured gas 
mass fraction with this radius, f gas (r25o o) = 0-129+p n}fi, is 
consistent with the value determined bv lAllen et al. {2008) 
from the earlier, shorter observation, but with significantly 
reduced statistical uncertainties. 



4. DISCUSSION 

We have presented new, deep Chandra observations of 
the high redshift X-ray cluster associated with 3C 186, 
a luminous comp act radio-loud quasa r at z = 1.067 
([Wills et al.lll992t ISchneider et alj|2007t) . The new obser- 
vations confirm the main re sults from our discovery paper 
(|Siemiginowska et al.ll2005l) . including the results on the 
global cluster temperature and central density profile. X- 
ray emission from the cluster is detected out to ~ 3 times 
larger distance from the quasar than was the case in the 
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Table 3 
Best Fit Model Parameters 



R° [arcsec] Range [a 



Total Counts 



Net Counts kT [keV] Norm c [le-3] CSTAT (dof=3550) n e [le-2 cm" 



3.375 


2.75-4.00 


603.0±24.6 


592.6±24.8 


3.111° 


91 

64 


37 1149+ 6 ' 0039 


2922.9 


4.875 


4.00- 5.75 


773.0±27.8 


751.8±28.2 


5.971J 


61 
25 


13.438311;™^ 


3064.0 


6.5 


5.75-7.25 


538.0±23.2 


513.9±23.7 


4.811^ 


61 
19 


7.4397l^; 1 1 9 1 f 6 


2976.1 


8.75 


7.25-10.25 


892.0±29.9 


827.1±30.9 


7.1ll^ 


43 
76 


9 R7O7 + 0-3014 
i.u 1 1 „ 2533 


3231.4 


12.75 


10.25-15.25 


1293.0±36.0 


1135.2±38.1 


7.771? 


93 
92 


1 nfi97 + °- 0996 
l.UOZ ( _ 0857 


3388.9 


18.25 


15.25-21.25 


1306.0±36.1 


1034.8±39.7 


6.95+? 


88 
44 


0.4663_ 0371 


3460.2 


25.625 


21.25-30.00 


1452.0±38.1 


896.4±44.8 


5.03l° 


6 7 

65 


n 97«1 +0-0153 
U.2781_„ 0156 


3533.8 



°- yi _0.45 

3.62l°;^ 

2 65+ 21 
1 59+ 00 

i ' u3 -0 .08 

1 00+ 05 

l.UU_ 04 
K + 0.03 
-0.03 
, +0.01 
-0.01 



0.66: 
0.51 4 



a The assumed annuli are circular with the mean radius listed in the R column and ranges in Range column; Total and net counts (in 0.5-7 keV 
range) summed within 4 observations in each region; c Normalization for APEC thermal model defined as Norm — - — ^ f n e riH dV with the 



abundance table set to Anders Grcvcssc (1989); d listed uncertainties are at 68% for one interesting parameter. 



Table 4 

Best Fit Parameters for the Simulated Quasar Data 



R a [arcsec] 



Norm 6 [l.c-7] 



3.375 
4.875 
6.5 
8.75 
12.75 
18.25 
25.625 



1 58+ 014 

i - os _0.17 
1 47+0.26 
-0.24 

± - '°-0.23 

1 46+ - 21 

, O- + 0.21 

i.oo_ 18 

1 - bS -0.35 



6.9^ 



3+0-9 
-0.8 

0.7 
- + 0.6 



2.5: 

3.1 
3.7 
1.5 



0.L 
+ 0.6 
-0.5 



7+0.7 
-0.6 
- + 0.5 
-0.4 
1 j+0.4 



a The assumed annuli are as in Table 3, e.g. circular with the centers listed in the R column; 
listed uncertainties arc at 68% for one interesting parameter. 
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Fig. 4. — Deprojected temperature, density and entropy profiles for the best fit model parameters shown in Table [3] 
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first, short observation. The quasar is located within the 
center of the diffuse X-ray emission and only slightly off- 
set (5.5±0.9 kpc) from the centroid of the cluster's X-ray 
emission. The high signal to noise data and a larger cluster 
area uncovered in the new observations allowed for more 
detailed analysis of the properties of the cluster. Below we 
discuss the main results of these new X-ray observations. 

4.f. Cooling core 

The 3C 186 X-ray cluster shows an elongated morphol- 
ogy that is detected out to r > 500 kpc. The cluster 
temperature profile has the characteristic shape of a cool- 
ing core cluster, with a sharp decline towards the center. 
The cluster is relatively cool in the outer (r > 200 kpc) re- 
gions with its temperature increasing slightly to a peak 
of 7.8±f£ kcV at 100-200 kpc and then declining to 
3.11±o 04 keV in the central regions (Fig. [§. The elec- 
tron density rises relatively smoothly from ~0.001 cm~ 3 
at r ~ 500kpc to ~ 0.1 cm~ 3 , as one moves inward from 
r=500kpc to the innermost resolved regions (Fig. [5]) . The 
sharp drop in entropy in the inner regions is also typical 
for cool core clusters. 

The cooling time profile is shown in Figure [6] We 
measure cooling times at radii of 50 kpc and 25 kpc of 
1.7±0.2 x 10 9 years and 7.5±2.6 x 10 8 years, respectively. 
We also measure nominal mass cooling rate (in the ab- 
sence of cooling) of 400 ± 190M Q year _1 within the central 
lOOkpc. The cluster surrounding 3C 186 clearly possesses 
a very strong cooling core. 

The measured core radius of ~ 28.6 ± 2.5 kpc is small in 
comparison to typical core radii of nearby clusters. How- 
ever, similarly small core radii have also been observed in 
other lower-redshift clusters with strongly cooling cores 
fg-g. lAllen etHI l2001t ISchmidt et all 120011: lAllen et all 
120021) . 

Obscrvationally X-ray clusters divide into two classes: 
cool core clusters and non-cool core clusters. This divi- 
sion is based on the p eak X- r ay su rface brightness and 
a central cooling time. iBurnsI (| 1990T) studied radio emis- 
sion of central galaxies (cD) in a sample of Abcll clusters 
and noticed that the cDs in cool core clusters were more 
likely than the ones in non-cool core clusters to be radio 
loud with a high radio power. More recent studies con- 
firm that the cool core clusters are more likely to host a 
cD with both radio emission and Ha emis si on lines (e.g. 
ISanderson et al.ll2009t iHaarsma et alj|2009t) . iMittal et all 
(|2009f ) reported that all strong cool core clusters harbor a 
central radio source, while 67% of the weak cool cores do, 
and 45% of the non-cool cores do. 

3C 186 cluster hosts a luminous quasar with broad lines 
and a compact radio structure fully contained within the 
host galaxy. 3C 186 is also a high rcdshift (z = 1.06) clus- 
ter with a str ong cooling core. This is contradictory to the 
suggestion bv lVikhlinin et al~l (|2007f ). based on a Chandra 
sample of X-ray clusters that there are no cool core clus- 
ters at z > 0.5. On the other hand more recent work 
([Santos et al.ll2008l : lAlshino et aLll2010D indicates that the 
fraction of weak to moderate cool cores remains the same 
at high redshift and only the fraction of strong cool cores 
drops significantly. 

Is 3C 186 unique? It is interesting to note that a fraction 
of AGN in clusters increases with redshift ([Martini et al.l 



l2009t ). so high rcdshift clusters are more likely to host an 
AGN. An X-ray emission associated with AGN can con- 
fuse detection of a cooling core in a cluster at high red- 
shift (|Branchesi et al.ll2007f ). Cool core clusters also have 
smaller cooling radii and their detection require high res- 
olution X-ray observations. Therefore, samples of X-ray 
clusters used in studies of cool core evolution with redshift 
may be missing clusters associated with a strong AGN. We 
note that H1821+643 cluster at z=0.3 also has a strong 
coolin g core that was observed by Chandra ( Russell ct al. 
l2010f ). Further studies of higher rcdshift clusters associ- 
ated with AGN are needed in order to understand the 
evolution of cluster cooling with redshift. 

4.2. Supermassive Black Hole Powering 3C 186 Quasar 

The black hole mass estimate for the 3C 186 quasar 
from meas urements of its broad emis sion lines is equal to 
3xl0 9 M Q (jSiemiginowska et al.ll2005l ). The corresponding 
Eddington luminosity is equal to L^dd — 4 x 10 47 erg s _1 . 
The quasar optical-UV lum inosity, based on the spectra l 
energy distribution given in ISiemiginowska et al.l (|2008al ). 
is equal to Lyy ~ 6 x 10 46 erg s~ 4 . U sing a bolometric 
correction that ranges between - 5 - 10 (jElvis et al.lll994h 
we estimate the 3C 186 bolometric luminosity to be of the 
order of Lb i ~ 10 47 erg s _1 with the required accretion 
rate of 0.25MEdd critical rate. 

The growth of this supermassive black hole might be 
closely related to the mass deposition from the cluster. If 
the cooling rate is ~ 470 M Q year -1 then only a small 
fraction, < 0.5%, of the cooling gas is needed to grow a 
10 9 M Q black hole within the cooling time of the cluster's 
core. However, the mechanism of transporting this gas to 
the close vicinity of a central black hole is unclear. 

4.3. Intermittent Radio Source 

3C 186 radio s ource belongs to a class of young CSS ra- 
dio so urces (see IO'Dealll998l for review). iMurgia et al.l 
(1999]) measured a synchrotron age of ~ 5 x 10 5 years 
for the entire radio structure. 3C 186 has a double radio 
morphology with a one-sided jet. Th e double radio source 
has a total length equal to 1.8" (see ISiemiginowska et al.l 
(2005) and references therein) corresponding to 15 kpc 
projected size. A deprojected size is at least 30 kpc for 
< 30 degrees angle to our line of sight (based on the one- 
sided VLA jet), but probably not much larger (e.g., 100 
kpc, or radius of 50 kpc, for ~ 9 degrees angle which is 
too small). Therefore we conclude that the 3C 186 radio 
source is contained within the host galaxy. 

Studies of compact radio sources suggest that they 
might have repe titive outbursts on short timescales of ~ 
10 3 - 10 5 years (iBaum et al" ll990HRevnolds fc Begelmar] 



Il997t lOwsianik et all 119981 : ICzernv et al.l I2009D . In "the 

case of the shortest timescales (~ 10 3 years) the radio 
source does not have enough energy to grow beyond the 
host galaxy and it starts to recollapse within the host glaxy 
ISM. If there have been previous outbursts of the radio ac- 
tivity in 3C 186 on the timescales longer than 3 x 10 4 years, 
the radio source would have been larger than the observed 
CSS structure. Our initial studies of the VLA radio data 
show a possible prese nce of an extended radio e mission 
on scales of 10 arcsec (jSiemiginowska et al.ll2008a|) . How- 
ever, this radio emission is seen at a very low significance 



10 



and more detailed analysis of new EVLA maps obtained 
recently have yet to to confirm that this structure is real. 

Most searches for X-ray clusters around radio-loud ac- 
tive galaxies have been focused on those with large-scale 
radio structures. Such radio structures are old, triggered 
a long time ago (> 10 7 years) and therefore have been in- 
teracting with the cluster environment for a long time. In 
nearby clusters, long term (~ 10 s years) intermittent radio 
activity of the central AGN is often imprinted into the X- 
ray morphology of the cluster in the form of bubbles, rip- 
ples or discon tinuities in the surface brigh tness indicative 
of shocks fsee lMcNamara fe Nulsen l 120071). GPS and CSS 
radio quasa rs are young (< 10 5 years iMurgia et a!1ll999t 
IO'Deafl998l for review) and have not developed large scale 
radio structures. These sources, if found in clusters, can 
potentially test the cluster heating process and the sig- 
nificance of the luminous quasar in the evolution of the 
cluster. 

4.4. Quasars and X-ray Clusters 

The majority of nearby clusters host a low power 
radio source with FRI radio morphology that have 
buoyantly rising bubbles filled with radio plasma. 
There are, however, a few examples of X-ray clus- 
ters associated with quasars or powerful radio galax- 
ies at lowe r redsh ifts (for example Cygn us A, 3C 295 
lAllen et al.1 (|2001l) . IRAS 09104+4109 llwasawa et all 
( 20011 ) . HS1821+643 lRussell et all (12010D ). We note that 
the FRII radio source would have pressure driven radio 
lobes and jets and its X-ray morphology may be different 
than the one seen in clusters with FRI radio sources. 

Cygnus A is a nearby FR II radio source em bedded 
in a bright X-ray cluster (|Carilli fe Barthell I1996D . The 
Chandra observations show a rich filamentary structure 
associated with the evolution of the radio source within 
the cluster and evi dence for heating of the cluster gas 
([Wilson et all I2000D . However, the Cy gnus A nucleus 
may not be in a luminous quasar pha se although its ra- 
dio power is high Lr ~ 10 45 erg s -1 (jYoung et al.l 120021 
ISteenbrugge et ai1l2008l ). The nucleus is hig hly absorbed, 
i.e. N h ~ 2 x 10 23 cm -2 as measured by lYoung et al.l 
(2002?), and correcting for the absorp tion gives the hard 
X-ray luminosity of 3.7xl0 44 erg s" 1 . lYoung et all (|2002t) 
estimated the optical luminosity to be consistent with 
Seyfert galaxies. 

Another case of a relatively bright X-ray cluster de- 
tected around a lower redshift z=0.32 2 luminous quasar is 
HS1821+643 (L h„i ~ 2 x 10 47 erg s" 1 iKolman et al.l fT99ll 
ICrawford et "all Il999t iRussell et all I2010D . This quasar 
has typical signatures of a quasar with broad lines and 
thermal emission i n optical-UV band. It hos ts a 300 kpc 
FRI radio source (|Blundell fc Rawlingsl 120011) that might 
be heating the cluster medium ( Russell et al.ll20Tfj|) . The 
Chandra observation indicates complex interactions be- 
tween the quasar and the cluster. However, the cluster 
properties are typical for the cool core cluster with a short 
central cooling time and the quasar does not appear to 
significantly impact the large scale cluster environment. 

4.5. Cluster Heating 



3C 186 is the first high redshift X-ray cluster known to 
host a luminous quasar and a compact radio source. The 
cluster X-ray morphology indicates that the cluster is well 
formed and has a cool core with a short central cooling 
time. The radio source can potentially supply the energy 
required to stabilize the cluster core against catastrophic 
cooli ng, as it expands into the c luster medium. 

In ISiemiginowska et al.l (120051) we estima ted the power 
of the radio jet using the lWillott et"aTI (|1999D (their Eq. 12) 
relation between radio luminosity and jet power defined 



Q 



3 x 10 38 L^ 5 l W, where L151 is in units of 



lO^W Hz" 1 sr" 1 . Assuming the 3C 186 151 MHz flux 
density of 15.59 Jy (lHales et al.lll993l ) which accounts for 
the total radio source emission (the radio core is absorbed 
at this frequency) we find L151 = 7.5 x 10 27 W Hz -1 sr -1 
and then from the above equation we obtain the jet power 
of L, et = 2.4 x 10 45 erg s -1 . There is significant scatter in 
the lWillott et al.l (|1999[ ) relation, therefore this is an order 
of magnitude estimate. 

The pressure in the radio lobes based on the radio 
flux density measurements and the equipartition ass ump- 
tion is ~ 10~ 8 erg cm -3 ([Sicmigin owska et al.ll2005l ) and 
it exceeds the thermal pressure of the clusteiQ ~ 4 x 
IQ-io erg cm -3 ( for kT = 31 keV and n = 0.O8 cm" 3 ). 

The overpressured radio source should drive a strong shock 
into the cluster medium and its expansion is not adiabatic. 
We can estimate a lower limit on the jet power from the 
equipartition measu rements. Using the radio lobes vol- 
ume of ~ 10 66 cm 3 (jSiemiginowska et al.l l2005h the mini- 
mum pressure gives a lower estimate of the instantaneous 
jet power, e.g. ~ 10 58 erg or ~ 6 x 10 44 erg s _1 for the 
age of the radio source of 5 x 10 5 years. 

We note that the bolometric luminosity of 3C 186 is 
equal to L hol ~ 10 47 erg s" 1 (see SecgTSJ) and the 3C 186 
radiative power exceeds the jet kinetic power by at least 
a factor of 10 . This is unusual for a cluster-center radio 
source and suggests that the so-called 'quasar mode' may 
be more important than the 'radio mode' for heating the 
3C 186 cluster. 

The process of transferring accretion energy into the 
cluster thermal energy is unclear. In the 'radio mode', the 
jet carries the energy from the black hole and deposits it 
into the cluster gas, e.g. via shocks. In the 'quasar mode', 
the radi ation s hould be a dominant carrier of the accretion 
energy. iKmel (120091) argues that for quasars with a black 
hole mass exceeding about 10 9 M Q , radiation can be very 
efficient in initiating strong outflows that cause radiative 
shocks and result in cluster heating. On the other hand, 
quasar radiation energy can also be transferred directly to 
the cluster gas via Compton scattering. 

We calculate the energy required to preve nt significan t 
cooling of the cluster core in 3C 186 following Eingl (|2009D . 
The mass of the gas within the cooling radius of 45 kpc 
assuming the central density of 0.08 cm -3 is equal to 
M C orc = 3.3 x 10 11 M Q . The amount of energy required 
to heat this gas is of the order of 1 keV per baryon, 
e.g. E heat ~ (lkcV/lGcV) M coro c 2 ~ 10- 6 M corc c 2 ~ 6_x 
10 59 erg. The observed quasar luminosity of ~ 10 47 erg s _1 
provides enough energy to heat the cluster core within 
about 2xl0 5 years if the heating process were 100 percent 
efficient. However, the cooling time of the core is much 



4 Note Sicmiginowska ct al. (2005) gives a factor of 10 lower value which is a mistake. 
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mass function is the most sensitive to the l inear growth 
of the fluctuations (jLinder fc Jenkins] 120031) . Thus the 
evolution of the number density of massive clusters that 
traces the growth of the density fluctuation can be used 
to constrain the dark energy equation of state parameter 
w (lHaiman et all 1200 lUVikhlinin et aUfeOOQl; iMantz et all 
120091) . The gas mass fraction of clusters can also be used to 
test the cosmological parameters (jAllen et al.ll2008l [2004) 
a nd provide indepen dent constraints on w. 

lAllen et al] (|2008l) measured f ga s(i"250o) f° r 42 clus- 
ters at < z < 1.1 and obtained an average value of 
f gas (r250o) = 0.1104 ± 0.0003. Their study included the 
previous, short observation of the 3C186 cluster (15ks 
good time) , which was the highest redshift cluster in their 
sample. The meas ured gas mass fraction for this target by 
lAllen et al.l d200l was / O a s (r 25 oo) = 0.1340 ±0.0777. Our 

new measurement of / sa s (72500) = 0- 129/^016 is consis- 
tent with the previous value, but improves the statistical 
uncertainties by a factor of ~ 4 — 5. Our new f gas mea- 
surement for the 3C 186 cluster is also consistent, within 
measurement errors, with the mean value determined at 
lower redshifts, arguing against any strong evolution of 
the /gas (72500) value for massive, relaxed clusters over the 
redshift range < z < 1.1. 

A relatively small number of X-ray clusters at z > 1 
have been observed so far (The BAX^ (|Sadat et al.ll2004f) 
data base lists 17 z > 1 clusters to date). New meth- 
ods to search for high-rcdshift X-ray clusters are providing 
new discoveries of massive clusters, including the recent 
discovery of the mos t distant X-ray cluster at z = 1.99 
(jAndreon et al.ll2009t ). The 3C 186 cluster is so far the 
only spectroscopically confirmed cooling core cluster at 
such high redshifts, and the only one known to host a 
luminous quasar at the center. The cluster global tem- 
perature, kT = 5.74 keV is not extreme when compared 
to the other massive clusters at z > 1, although the X- 
ray luminosity, L .5- 2 kcV = 6.4 x 10 44 erg s _1 , exceeds 
that of the next most luminous cluster d etected so far, 
XMMU J2235. 3-2557 dRosati et al. 1 I2009D at z = 1.39. 
The measured 3C 186 cluster mctallicity lie s slightly above 
the m ean trend with redshift discussed by iBalestra et al.l 
( 2OO7]) . More detailed evolutionary studies should soon be- 
come possible with improved high redshift cluster samples 
from Chandra and XMM-Newton. 

5. SUMMARY AND CONCLUSIONS 

We have presented results from a deep (200ks) Chan- 
dra image of the hot (kT = 5.6 ± 0.3kcV), X-ray luminous 
galaxy cluster surrounding the powerful quasar 3C186 at 
a redshift z = 1.067. The spatially resolved temperature 
profile, entropy, density and cooling time profiles all con- 
firm 3C 186 as a cooling core cluster. The measured gas 
mass fraction of f gas = 0.129tQ'o];g at ^2500 is consistent 
with measurements for lower redshift systems. This argues 
against strong evolution in f gas (z) at f25oo for massive, re- 
laxed systems. Cooling gas in the cluster core can in prin- 
ciple support the growth of a supermassive black hole and 
power the luminous quasar. The radiative power of the 
quasar exceeds the kinematic power, suggesting that ra- 

5 We note that the UV photons will also cool the cluster gas, however the total energy contained in the quasar spectrum above > IkeV exceeds 
the energy in the UV band. The heating-cooling balance of the cluster gas has to be modeled properly to understand the heating efficiency. 

e http://bax.ast.obs-mip.fr/ 



longer, ~ 7 x 10 8 years, and only a very small fraction of 
the observed luminosity is needed to support the cluster 
heating within that time. 

We estimate the efficiency of transferring the radia- 
tion energy into the cluster. The opacity of the clus- 
ter gas in the core to Compton scattering is equal to 
r = ii e rQ(7T ~ 0.008, assuming the average particle den- 
sity of n e = 0.075 cm~ 3 , r = 45 kpc core radius, and 
<tt is the Thomson cross section. This means that about 
0.8% of the quasar photons will interact with cluster gas 
and may heat up the cluster. We note also that only the 
photons with energies higher than ~ 1 keV will be able to 
hcatQ the cluster gas to the required temperatures which 
reduces the amount of available photons by a factor of 10. 
We also assumed spherical symmetry and if we include a 
covering factor of 30% that accounts for cold gas in the 
central regions (e.g. torus), the total number of photons 
decreases to about 0.3%. This results in the available in- 
stantenous power of 0.003 x 10 46 erg s _1 which is equivalent 
to a total energy of ~ 7 x 10 59 erg available to balance the 
cluster cooling within 7x 10 8 years. This is within an order 
of magnitude of the required energy to balance the cluster 
cooling. However, this also requires that the quasar is pow- 
ered by accreting at 0.25Msdd rate, (e.g. ~ 18M© year -1 
assuming standard 10% accretion efficiency), continuously 
during the cooling time of the cluster, which results in the 
final required amount of fuel exceeding by a factor of 4 the 
current black hole mass. Although such an amount of fuel 
is only a fraction of the available mass in the cluster core, 
it is unclear how it can be funneled down to the center 
w ithin close v icinity of a black hole. 

iKind (|2009D suggested that the radiation pressure in the 
luminous quasars will generate a wind, so the radiation en- 
ergy is transferred to the kinetic energy of the wind and 
then to the cluster gas via radiative shocks. This mecha- 
nism also requires a large amount of fuel supply. 

We should comment here that the young CSS radio 
sources, such as 3C 186, might be intermittent on short 
timescalcs ~ 10 5 year s if the accretion r ates are close to 
the Eddington value. (ICzernv et al.l |2009f). The source ex- 
periences transitions between (1) a high super-Eddington 
state characterize by a luminous accretion disk and a pow- 
erful jet; and (2) a quiescent state with the sub-Eddington 
luminosity of the disk and no jet. In this scenario the 
super-Eddington state should last longer than ~ 10 3 years 
for a radio source to grow beyond the host galaxy, as the 
expansion of the radio source within the host galaxy takes 
more than ~ 10 6 years. During that initial time the lumi- 
nous quasar can provide enough radiative power to heat 
up a small cluster core. At later times when the radio 
source grows beyond the host galaxy and expands within 
the cluster environment the mechanical energy is directly 
used to prevent cluster cooling. 

4.6. Clusters at High Redshift and Cosmology 

X-ray observations of galaxy clusters at high redshift 
can provide strong co nstraints on cosmological parameters 
(see IVikhlininl (|2010| ) for a recent review). In hierarchi- 
cal models of cluster formation, the high-mass end of the 
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diativc heating my be important at intermittent intervals 
in cluster cores. 
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APPENDIX 

Al: Modeling Background Spectra 

We used Sherpa to fit simultaneously source and background data in the spectral analysis of the cluster. A "blank- 
sky" background file acis7sD2005-09-01bkgrnd_ctiN0001 . f its provided by the Chandra X-ray Center was used to 
define a background model. First we filtered and reprojected X-ray events contained in this background file follow- 
ing the CIAO Thread The ACIS "Blank-Sky" Background File$H Next we extracted a spectrum from a box region 
(rotbox (4385. 5, 41 12. 5, 980, 990,0), the first two values show a center of the box, the next ones the size of the box in 
ACIS pixels, and a rotation angle) and corresponding instrument response files using specextract tool in CIAO. We fit 
this spectrum in Sherpa using a combination of an 8th-order polynomial and 9 gaussian lines getting the best-fit model 
with the statistics equal to cstat=540.8 (443 d.o.f). Figure [7] shows the resulting fit and residuals. 

We next checked how well this empirical model fits the background data in the observations. We applied the background 
model to the background spectra in each obsid and fit only the model normalization. The resulting fit statistics is equal 
to cstat=1900.4 (1772 d.o.f.) for the simultaneous background fit to all four observations. Figure [8] shows the fit result 
for the obsid 9408 with the highest background counts in the spectrum. 

In the final simultaneous source and background fitting of the cluster spectra the background normalization was varied 
and appropriate background model predicted counts were included in the total model predicted counts. An additional 
constant scaling accounts for a difference in exposure times and areas of the background and source regions and it is 
automatically applied by Sherpa during the fit. 



http: / / cxc.harvard.edu / ciao / threads / acisbackground / index.py.html 



14 



Blank-Sky Background 




1 



Energy (keV) 



Fig. 7. — Top panel: Blank-sky background data fit with the complex model: polynomial + 9 gauss lines. Lower panel: 
Residuals. 
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Fig. 8. — The background spectrum for obsid 9408 fit with the same shape background model as the blank-sky data shown in 
Fig. Top panel shows the data marked with points and a model drawn with a red solid line. The bottom panel shows the 
residuals. The data were grouped by 15 counts per bin for the visualization only. 



